close all;
clear all;
clc;

%define grid
Sx=1;  %physical size along x
Sy=1;  %physical size along y
Nx=10; %number of cells along x
Ny=10; %number of cells along x


% grid arrays
dx = Sx/Nx;
xa=[0:Nx-1]*dx;
xa = xa - mean(xa);

dy = Sy/Ny;
ya=[0:Ny-1]*dy;
ya = ya - mean(ya);


% create circule
r = 0.4;
[Y,X] = meshgrid(ya,xa);

A = (X.^2 + Y.^2) <= r^2


